Influence of the Third Dimension of Quasi-Two-Dimensional Cuprate 
Superconductors on Angle-Resolved Photoemission Spectra 
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Angle-resolved photoemission spectroscopy (ARPES) presents significant simplications in analyz- 
ing strictly two-dimensional (2D) materials, but even the most anisotropic physical systems display 
some residual three-dimensionality. Here we demonstrate how this third dimension manifests itself 
in ARPES spectra of quasi-2D materials by considering the example of the cuprate Bi2Sr2CaCu2 08 
(Bi2212). The intercell, interlayer hopping, which is responsible for fcz-dispersion of the bands, is 
found to induce an irreducible broadening to the ARPES lineshapes with a characteristic dependence 
on the in-plane momentum fcy . Our study suggests that ARPES lineshapes can provide a direct 
spectroscopic window for establishing the existence of coherent c-axis conductivity in a material 
via the detection of this new broadening mechanism, and bears on the understanding of 2D to 3D 
crossover and pseudogap and stripe physics in novel materials through ARPES experiments. 

PACS numbers: 79.60.-1, 71.18.+y, 74.72. Hs 



Angle-resolved photoemission spectroscopy (ARPES) 
has been applied extensively in the recent years for in- 
vestigating the electronic and quasiparticle properties of 
high-temperature superconductorsi. Much of the exist- 
ing ARPES work on the cuprates and other quasi-2D 
materials implicitly assumes perfect two-dimensionality, 
ignoring the effects of dispersion in the third dimension. 
As energy resolutions of the order of a few meV's have 
now become possible in the state-of-the-art ARPES in- 
strumentation, it is natural to ask the question: How 
will /cz-dispersion play out in the ARPES spectra of the 
quasi-2D materials? An obvious answer is that spec- 
tral peaks will undergo shifts with photon energy -albeit 
small- much like in the 3D case^*^ Our analysis however 
reveals a very different scenario in that insofar as ARPES 
response to fcz-dispersion is concerned, quasi-2D systems 
differ fundamentally from their 3D counterparts. We find 
that in the presence of typical final state dampings of the 
order of eV's, initial state dispersions over much smaller 
energy scales will not appear as energy shifts with hv, 
but will instead induce an irreducible linewidth in the 
spectral peaks with a characteristic k\\ dependence. 

Our results bear on a variety of important issues in 
cuprate physics and give insight into a number of puz- 
zling features of ARPES in the cuprates. Since our new 
broadening mechanism does not have its origin in a scat- 
tering process, it explains why the lifetimes derived from 
ARPES spectra are generally found to be shorter than 
those obtained from other experiments. Moreover, Fermi 
surface (FS) maps obtained via ARPES often display 
broad patches of spectral weight rather than well-defined 
FS imprints, particularly near the antinodal point; such 
patches occur naturally in our calculations. Indeed, much 
of the pseudogap phenomena - in particular the lack of 
well-defined quasiparticles near (tt, 0) - will need to be 
reevaluated in light of the present findings. In this vein, 
these results also impact the analysis of stripe physics 
through ARPES experiments. Finally, detection of the 



related linewidths in the ARPES spectra offers a unique 
spectroscopic window for establishing the existence of co- 
herent c-axis conductivity and intercell coupling in a sys- 
tem; this important intrinsic property is not accessible 
directly through other techniques,"'^ 

Specifically, we focus in this article on the tetrago- 
nal body-centered Bi2Sr2CaCu208 (Bi2212) compound, 
which has been a workhorse of ARPES studies. Bi2212 
is a nearly 2D material with two Cu02-layers in the 
primitive unit cell. The intracell interaction between 
the two Cu02-layers, spaced a relatively short distance 
of ~ 3.2 A apart, results in the well-known bilayer 
splittingiiSiL2i2ii£. On the other hand, the intercell cou- 
pling between the bilayer slabs in different unit cells is 
expected to be smaller due to the larger intercell Cu-Cu 
distance of ~ 12 A resulting in weak but non-vanishing 
fc^-dispersion. The intercell coupling will be even more 
pronounced in other high-Tc's since Bi2212 presents one 
of the longest c-axes in the cuprate family. 

Concerning computational details, the band struc- 
ture for Bi2212 was obtained within the local-density- 
approximation (LDA) by using the well-established 
Green's function methodologjiii. The crystal potential 
used is the same as that employed in our previous stud- 
ies of Bi2212 and involves 30 atoms per conventional unit 
cell^i^i^, yielding good agreement with the experimen- 
tally observed Fermi surface (FS)iaiiiiii^. ARPES 
intensities have been computed within the one-step pho- 
toemission formalism; see, Refs. 1^ andlvl for details. 

Figure ^a), which considers the familiar antibonding 
(A) and bonding (B) bands in Bi2212, shows that the 
k|| -dispersion depends strongly on kz as the associated 
bands wander over the two sets of shaded areas. The k^- 
dispersion (i.e. the vertical width of the shaded areas) 
displays a striking dependence on fcy and nearly vanishes 
at the antinodal point fcy = (tt/a, 0). A clear bilayer 
splitting between the A and B bands is seen at all de- 
values, except at fcii = fcit = 0.2(27r/a) for A:^ = (solid 
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FIG. 1; (a): Calculated first-principles ky -dispersion in 
Bi2212 along [100]-direction for the B and A bands at three 
difi'erent values (in units of 2it/c): kz — (solid lines); 
kz = 0.5 (dashed), and kz = 1 (dash-dotted). Shading de- 
notes the regions over which the bands wander as a function 
of kz- (b): fez-dispersion of the B band at five difi'erent fey 
values ranging from to vr/a. (c): Same as (a), except that 
these results are based on a tight binding formalism. 



lines), where A and B bands touch. This level crossing 
leads to the anomalous fcz-dispersion shown in Fig. 1(b) 
for the B band: The shape changes from being described 
approximately by the form sin^(A:zc/4) for fcy = 0, to 
|sin(fczc/4)| at fcj|, and finally back to sin^(fczc/4) for 
/c||a/(27r) = 0.5. The antibonding band behaves simi- 
larly (not shown for brevity). The complex behavior of 
the bilayer splitting and the fc^-dispersion shown in the 
first-principles computations of Figs. 1(a) and 1(b) can 
be modeled reasonably well by a relatively simple tight- 
binding (TB) Hamiltonian as seen from Fig. 1(c). We 
return below to discuss this point further. 

We next turn our attention to how fcz-dispersion will 
manifest itself in ARPES spectra. The key insight in 
this regard is provided by the simulations^^ of Fig. 2, 
which consider the evolution of the lineshape of the en- 
ergy distribution curves (EDCs) as the final state broad- 
ening introduced via the imaginary part of the final state 
self-energy^^, Sy, is increased from a very small value 
in (a) to a realistic value in (c). In (a), the position 
of the spectral peak undergoes the familiar shift as hv 
varies. This shift results from the fact that in the photo- 
excitation process /c|| must remain unchanged in trans- 
mission of the electron across the surface and the initial 
and final states can connect only at a specific value of 
kz in order to conserve energy. Note that the total shift 
in the peak position in the EDCs of (a) gives the size of 
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FIG. 2: Simulated ARPES hneshapes (EDCs) in Bi2212 for 
a series of photon energies (hf — 25 — 25.7 eV) at a fixed 
fe|| — (0.34, 0.09)27r/a-point using three different values of the 
final state broadening given by the indicated imaginary parts 
of the self-energy, E". In order to highlight the infiuence of 
fez-dispersion, the initial state broadening is chosen to be very 
smaU, E" = 0.2 meViii 



the fcz-dispcrsion of the initial state. This should not be 
confused with the change in hv needed to probe such a 
band. The change in hiy is controlled by the final state 
dispersion, which is generally much larger than that of 
the initial state. 

Continuing to the intermediate case of =0.1 eV in 
Fig. 2(b), the shift in the peak position from the bottom 
to the top of the initial state band is once again evident, 
but the lineshapes are quite different, even though the 
initial state damping S" in (b) is identical to that in 
(a). It is striking that some spectral intensity appears 
in (b) at all energies encompassed by the initial state 
band at every hv. This remarkable effect comes about 
because the energy uncertainty permitted by the width 
of the final state allows the photoelectron to couple with 
initial states off-thc-energy-shell. The lineshape thus de- 
velops a new component with an ft,i/-independent width 
equal to the initial state bandwidth, which rides on top 
of the energy conserving peak in (a). Notice also how 
the changes in the peak position could allow observing 
dilferent FSs as one maps different values of k^- Finally, 
for the realistic final state width of = 1 eV in (c), the 
initial state bandwidth component dominates the line- 
shape. The lineshape of the EDC curve is now virtually 
/ii'-independent. Despite the large final-state broaden- 
ing, the energy spread of the EDC remains equal to the 
initial state bandwidth in kz-direction because outside of 
this interval, there are no initial state electrons capable 
of absorbing the photon. 

Fig. 3 elaborates on these points by considering some 
of the results of Fig. 2 over a much broader range of pho- 
ton energies. Fig. 3(a) presents the spectra of Fig. 2(b) 
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FIG. 3: (Color) (a) Simulated ARPES intensities of Fig. 2(b) 
are depicted over a wide photon energy range of 19 — 27 eV 
for E^' = 0.1 eV. Reds denote highs, (b) Final state band 
structure as a function of fcz. (c) Same as (a) except that 
E" = 1 eV, which corresponds to the case of Fig. 2(c). 



for = 0.1 eV for /ii/ w 19 - 27 eV. Different colored 
bands here represent the same initial states being excited 
to different final states. Comparing Figs. 3(a) and 3(b), 
we see that the peaks of Fig. 3(a) more or less follow the 
final state band structure of Fig. 3(b) [The agreement is 
not expected to be perfect of course, since the initial state 
dispersion is not entirely negligible]. In Fig. 3(c), when a 
more realistic final state width of = 1 eV is assumed, 
much of the structure is lost. The large variations in col- 
ors (or intensities) observed in Figs. 3(a) and 3(c) are the 
consequence of the well-known fc|| and hi^ -dependency of 
the ARPES matrix element 

The preceding discussion of Figs. 2 and 3 makes it ob- 
vious that in a quasi-2D system, fcz-dispersion leads to an 
irreducible linewidth in the ARPES spectra, which can- 
not be resolved by changing photon frequency. This effect 
will also appear in the FS maps observed by measuring 
emission from Ep in the {kx, A:y)-plane. Fig. 4 illustrates 
how this plays out. The standard bilayer split FS in the 
kz = plane, given in Fig. 4(a), is the FS usually thought 
to be measured in ARPES. In Fig. 4(b) the fuU 3D FS 
has been projected on to the (OOl)-plane by collecting 
individual FS cuts corresponding to different k^ values. 
The fc^-dispersion is now seen to introduce an effective 
"broadening" to the FS imprint. In a perfect 2D-system, 
the maps in (b) and (a) would be identical. The effect of 
using a finite energy window of ±30 meV around Ep in 
the computations is shown in Figs. 4(c) and 4(d). These 
simulations indicate the influence of a finite experimen- 
tal resolution on the results, Note that even with the 
window of ±30 meV in (d), the broadening effect of kz- 
dispersion is not washed out. In fact, the broadening is 
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FIG. 4: Different imprints of the FS of Bi2212. (a): Section 
of the FS at fcz = 0. (b): Projection of the FS on to the 
(OOl)-plane, where the "broadening" of the FS reflects the 
effect of fez-dispersion, (c) and (d) correspond to (a) and (b), 
respectively, except that the effect of adding a finite energy 
window of ±30 meV is included in the computations. 

somewhat enhanced, especially near the antinodal point 
due to the contribution of the flat bands related to the 
van Hove singularity (VHS). 

It is also interesting to consider changes in the 
linewidth in momentum, Ak, as one moves away from 
the Fermi level. In general, our simulations indicate that 
Afc increases with increasing binding energy (BE), due 
mainly to the flattening of bands and the concomitant 
reduction in the band velocity.^ In any event, the broad- 
ening as a function of BE is neither simply quadratic nor 
exponential. Furthermore, considering the the effect of 
varying S" at a general fc|| -point, we have found that as 
the initial state damping due to intrinsic scattering mech- 
anisms (simulated via the value of E") decreases, the 
linewidth AE becomes increasingly dominated by the ir- 
reducible width associated with /c^-dispersion. Along the 
antinodal direction this dispersion is negligible, and AE 
and S" are related linearlyi^i 

We return now to comment briefly on the bands of 
Figs. 1(a) and 1(b). Insight into how the complex fcy- and 
/cz-dependencies of these bands reflect intercell as well 
as intracell hopping effects can be gained by modeling 
these bands within the TB framework. In the absence 
of intercell coupling, the conventional bilayer splitting 
possesses the form = tz{cx — , with Ci = cos{kia), 
i ~ x,y. The intercell coupling in the cuprates may 
be included in the one-band model with dispersion — 

ek = ~2t{Ca; + Cy) - WCxCy - ( fc 1 1 , ) [ (Cj; " Cy)^/4± flo], 

where Sz — sin{kzc/4:) and ag = 0.6 corrects for the 
finite bilayer splitting found at /cy — 0. The form 
of Tz depends on the particular cuprate considered. 
While Ref. lists seven inequivalent interlayer hopping 
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parameters, we find that we can describe the dispersion 
reasonably by including only the dominant contribution 
associated with hopping between C u 4s levels by intro- 
ducing Tz — ±-\/ {tz — t'zY + '^izt'zsl, where plus (minus) 
sign refers to the bonding (antibonding) solution, is 
a constant associated with intracell interlayer hoppings, 
and t'z = At'zQCOs{kxa/2)cos{kya/2) accounts for intercell 
hopping. The extra angular dependence in this term 
arises because for intercell hopping the two Cu02-planes 
are offset, so that one Cu atom sits above an empty 
cell. Note that for = one obtains the simple bilayer 
sphtting tbi with no fc^-dispersion. TB bands of Fig. ^c) 
assume (in eV): t = 0.42, t' = -0.12, tz = 0.125, and 
t^g = 0.04. A comparison of the first principles and TB 
bands in Figs.^a) and^c), respectively, shows that our 
TB model reproduces the bilayer splitting, the collapse 
of the fc^-dispersion around (tt/o, 0), and the anomalous 
dispersions associated with level crossing to a reasonably 
good degree, some minor discrepancies notwithstanding. 

We emphasize that t'^,, which provides intercell cou- 
pling, is essential for obtaining coherent c-axis conduc- 
tivity, t'z thus controls the intrinsic resistivity anisotropy 
and gives a measure for discriminating between coher- 
ent and incoherent c-axis hopping^ The lifetimes in the 
cuprates extracted from ARPES are considerably smaller 
than optical lifetimes or those deduced from transport 
or tunneling measurements^^. A part of this difference 
may be explained to be the result of an unresolved bi- 
layer splittingSi. The remaining anomalous broadening 
is often interpreted^^ in terms of small-angle scattering, 
which does not contribute to transport, but its origin re- 
mains obscureSi. The present analysis indicates that part 
of the broadening of ARPES lines has its origin in the k^- 
dispersion, unrelated to any scattering mechanism. 

Our predicted effects of A:2-dispersion on the 
linewidths in Bi2212 should be resolvable with the 
power of currently available high resolution ARPES 
instrumentation^ The size of our novel line broadening 



mechanism will be larger in materials with greater 3D 
character (e.g. YBCO, LSCO, NCCO, ruthnates, man- 
ganites, etc.), although we expect the results of Figs. 1-4 
to provide at least a qualitative handle in quasi-2D mate- 
rials more generally. Notably, coherent 3D coupling has 
recently been reported in a Tl cuprate^^. Additionally, a 
similar fcz-induced broadening may be expected to have 
an impact on other experimental probes, such as Raman 
scattering2£. 

In conclusion, we have demonstrated that residual k^- 
dispersion in quasi-2D materials will induce an irreducible 
linewidth in ARPES spectra. This intrinsic linewidth of- 
fers a new spectroscopic window for understanding in- 
as well as out-of-plane scattering mechanisms and the 
nature of the 2D to 3D crossover in the cuprates. This 
highly anisotropic line broadening mechanism, which has 
not been recognized previously and is unrelated to 2D 
physics, indicates that the existing analysis of stripe and 
pseudogap physics based on ARPES spectra should be 
reexamined. Our study shows how ARPES can be ex- 
tended to unravel the hidden third dimension of energy 
bands and Fermi surfaces in quasi-2D systems with wide- 
ranging consequences for understanding the nature of 
electronic states in many novel materials. 
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